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The present invention provides a practical method for UAVs 
to hike advantage of thermals in a maimer similar to piloted 
aircrafts and soaring birds. In general, the invention is a 
method for a UAV to autonomously locate a thermal and be 
guided to the thermal to greatly improve range and endurance 
of the aircraft. 
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Figure 1. Top level block diagram of the present invention. 
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Figure 2. Block diagram of the total energy estimation (block 10). 
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wup 

Figure 3. Block diagram of the thermal identification (block 20). 



Figure 4. Comparison of the queue, q 5 , before the drift correction to the queue, q c , after 
the drift correction. Drift was calculated using P up/ and P up 2 - 
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Figure 5. Predicted updraft velocity distribution from equation 17 for w up =2m/s and 
r up = 45m. 



Figure 6. Influence vector, H, for a queue size, N q , of 45. The influence vector is used to 
de-emphasize older data points. 
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Figure 7. Block diagram of soaring controller (block 40). 
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GUIDANCE AND CONTROL FOR AN 
AUTONOMOUS SOARING UAV 

STATEMENT OF GOVERNMENT INTEREST 

5 

The invention described herein may be manufactured and 
used by or for the Government of the United States of 
America for governmental purposes without payment of any 
royalties thereon or therefor. 

10 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates generally autonomous soar- 
ing for uninhabited air vehicles (UAVs), particularly to 15 
improved endurance for autonomous soaring UAVs, and 
more particularly to using thermals or convective thermals to 
improve the performance of UAVs. 

2. Description of the Related Art 

Glider pilots have been using buoyant plumes of air found 20 
in the lower atmosphere, normally referred to as thermals, in 
order to extend the distance that a glider may travel for many 
years. 

Several systems have been developed to assist pilots of 
gliders or low-powered aircraft to locate such thermals in 25 
order to provide improved range, endurance, or cross-country 
speed. For example, U.S. Pat. No. 6,012,675 describes a 
system that monitors air humidity to locate “moist” thermals 
and U.S. patent application 2004/012907 1 describes a system 
that employs a database of places on the ground that are likely 30 
to generate thermals in conjunction with current weather 
projections to attempt to estimate positions of thermals for 
pilots. 

Over the past several years, it has also been suggested that 
UAVs could benefit from using thermals in the same manner 35 
as piloted aircraft to improve range, endurance, and speed. 

For instance, Wharington, et al. in “Control of High Endur- 
ance Unmanned Air Vehicle” proposes the concept of using 
thermals to improve the performance of UAVs and discloses 
a general control scheme that could be applied for the UAV to 40 
take advantage of thermals. However, the paper states that a 
superior control scheme would be necessary for practical 
application of the concept. Further, Michael J. Allen, the 
inventor of the present invention, in “Autonomous Soaring for 
Improved Endurance of a Small Uninhabited Air Vehicle” 45 
provides mathematical evidence that UAVs may dramatically 
improve range and endurance when employing thermals. 

However, until the present invention, no known practical 
control and guidance system has been developed for UAVs to 
take advantage of the benefits of thermals. Therefore, it is 50 
desired to provide a guidance and control system and method 
for UAVs to locate and use thermals to provide improved 
range and endurance. 

SUMMARY OF THE INVENTION 55 

The invention proposed herein comprises a system and 
method of guidance and control for autonomous UAVs that 
allow the UAVs to locate and use thermals in order to increase 
the range and endurance of the UAVs. 60 

Accordingly, it is an object of this invention to provide a 
system and method for an autonomous UAV to independently 
locate thermals in the lower atmosphere. 

It is a further object of this invention to provide a system 
and method for guiding an autonomous UAV to the “center” 65 
or “strongest” portion of a located thermal in the lower atmo- 
sphere. 


2 

This invention meets these and other objectives related to 
improving guidance and control for an autonomous soaring 
UAV by providing a method for locating and using thermals 
to improve endurance of UAVs. The method includes the step 
of estimating energy rate and energy acceleration from total 
energy of the air vehicle. This is done by using values for 
static pressure around the UAV, true air speed of the UAV, and 
command to the motor. The method also includes the step of 
identifying the radius, strength and position of a thermal. This 
is accomplished using the energy rate and location of the 
UAV. The UAV position error, velocity error, and steady-state 
turn rate are then obtained using the location of the UAV, and 
the thermal radius, strength, and position. The method also 
includes identifying the soaring turn rate command for the 
UAV using the UAV position error, velocity error, steady- 
state turn rate, and energy acceleration. Also, a final turn rate 
or bank angle is identified for the UAV through mode switch- 
ing using the soaring turn rate, and a waypoint tracking turn 
rate. Finally, the UAV is controlled to guide it to locate and 
stay within the thermal. 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the drawings, which are not necessarily to scale, like or 
corresponding parts are denoted by like or corresponding 
reference numerals. 

FIG. 1 depicts a top level block diagram of the present 
invention. 

FIG. 2 depicts a block diagram of the total energy estima- 
tion shown in FIG. 1. 

FIG. 3 depicts a block diagram of the thermal identification 
shown in FIG. 1. 

FIG. 4 depicts a graph comparing the queue before and 
after drift correction is applied as set forth in FIG. 3. 

FIG. 5 depicts a graph showing predicted thermal velocity 
distribution from equation 1 7 in the specification. 

FIG. 6 depicts an influence vector used to deemphasize 
older data points in the queue shown in FIG. 3. 

FIG. 7 depicts a block diagram showing an embodiment of 
the controller of the present invention. 

FIG. 8 depicts a block diagram showing the mode switch- 
ing logic employed in the present invention. 

DESCRIPTION OF PREFERRED 
EMBODIMENTS 

The following nomenclature will be employed herein to 
describe the present invention: 
a equatorial radius of the earth, m 
D estimated thermal drift velocity, m/s 
err^, position error, m 
err sumsq sum squared error, m 2 /s 2 

eir su sum squared error for perturbed thermal radius, 
m 2 /s 2 

err,, velocity error, m/s 
err w weighted error, m/s 

err w weighted error for perturbed thermal radius, m/s 
E energy rate, Nm 

E motor to t a l energy rate due to the motor, N-rn/s 

E energy acceleration, nr/s 2 

E„ normalized energy acceleration, m/s 2 

g gravitational constant, m/s 2 

H influence vector 

h est estimated aircraft altitude, m 

i queue index 

k current measurement index 
k : energy acceleration gain, unitless 
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k 2 position error gain, unitless 
k 3 velocity error gain, unitless 
L learning rate, 1/s 2 
Lat latitude of aircraft, rad 

Lat 0 latitude of X Y reference frame, rad 5 

Lon longitude of aircraft, rad 
Lon 0 longitude of X Y reference frame, rad 
M aircraft mass, kg 
N number of rows in the queue. 

^static static pressure, Pa to 

P estimated thermal position, m 

P upDl thermal position calculated from the first 20 entries in 
the queue, m 

P upD thermal position calculated from the last 20 entries in 
the queue, m 15 

P^, aircraft X position, m 
P r aircraft Y position, m 
q : queue of aircraft position and energy rate 
q^ queue of aircraft position and shifted energy rate 
q 2 queue of drift corrected aircraft position and energy rate 20 
r estimated thermal radius, m 
r up perturbed thermal radius, m 
s laplace transform variable 

S distance between the aircraft and the thermal center, nt 

soarMode mode logic logical 25 

T thrust from the motor, N 

T max maximum thrust from the motor, N 

throt^! one-second delayed throttle command, fraction 

throt rrim trim throttle needed for level flight, fraction 

V airspeed, m/s 

V motor aircraft velocity rate due to the motor, no/s 2 
w up estimated thermal velocity, m/s 
w predicted thermal velocity, m/s 

'bpred predicted thermal velocity for perturbed thermal 35 
radius, m/s 
8,. stepsize, m 

A,, error gradient for perturbed thermal radius, m/s 2 
turn rate command, deg/s 

steady state turn rate, deg/s 40 

One relatively unexplored way to improve the range, dura- 
tion, or cross-country speed of an autonomous aircraft is to 
use buoyant plumes of air found in the lower atmosphere 
called thermals or convective thermals . Thermals occur when 
the air near the ground becomes less dense than the surround- 45 
ing air because of heating or humidity changes at the Earth’s 
surface. Piloted sailplanes rely solely on this free energy to 
stay aloft for hours at a time and to make cross-country flights. 
Soaring birds, such as hawks and vultures, have been 
observed to circle for hours without flapping their wings, and 50 
Frigatebirds are known to soar continuously day and night 
using thermals. 

Many UAVs have similar sizes, wing loadings, and mission 
profiles to soaring birds and sailplanes. Mission profiles that 
could allow small UAVs to take advantage of thermals include 55 
remote sensing, surveillance, atmospheric research, commu- 
nications, force protection, and Earth science. A study using 
a simple UAV simulation with thermals calculated from mea- 
sured surface and balloon data found that a 2-hour nominal- 
endurance UAV can gain up to 12-hours of flight time using 60 
thermals. 

The present invention provides a practical method for 
UAVs to take advantage of thermals in a manner similar to 
piloted aircrafts and soaring birds. In general, the invention 
comprises a method for a UAV to locate a thermal and autono- 65 
mously be guided to the thermal to greatly improve the range 
and endurance of the aircraft. 


The top-level autonomous soaring guidance and control 
method is shown in FIG. 1. Static and total pressure sensor 
input (to obtain true air speed) along with the aircraft throttle 
command were used to determine the rate (first derivative) 
and acceleration (second derivative) of aircraft total energy in 
the total energy estimation block 10. These energy rate and 
energy acceleration values are used as inputs to the mode 
switching block 50 as described further below. 

The energy rate from the total energy estimation block 10 
along with the position of the aircraft (latitude and longitude) 
are then used to identify a thermal that has been encountered 
by the aircraft 20. The thermal radius, strength, and position 
are calculated in block 20, using these input values. 

In the circle guidance block 30, the thermal radius, 
strength, and position from block 20 are used to obtain values 
for aircraft position error and aircraft velocity error in relation 
to the identified thermal position. The UAV steady-state turn 
rate is also obtained in this step. 

The aircraft position error, aircraft velocity error, aircraft 
steady-state mm rate acquired in guidance block 30, and the 
energy acceleration obtained from the total eneigy estimation 
block 10 are used to provide a soaring mm rate command in 
the controller block 40. The soaring mm rate command is the 
turn rate that is given to the UAV by the controller to lead the 
aircraft within the thermal to maximize climb rate while the 
UAV is in soaring mode. 

Finally, the energy rate, energy acceleration, waypoint 
tracking turn command, and soaring mm rate command are 
input into the mode switching block 50. In the mode switch- 
ing block 50, mode logic is used to determine when the UAV 
should be told to switch between searching flight (looking for 
a thermal) and soaring flight (using a thermal). Energy rate 
and energy acceleration are used to determine when the air- 
craft should be in searching mode and follow the waypoint 
tracking turn command or if the aircraft should be in soaring 
mode and follow the soaring rate turn command. In this block, 
when the amount of energy obtained from a thermal decreases 
to a certain level, the UAV is commanded to switch from 
soaring flight to searching flight. 

FIG. 2 depicts a more detailed explanation of the total 
energy estimation block 10 of FIG. 1. Aircraft total energy 
(non-engine energy) 200 is first calculated from airspeed, V, 
and static pressure, P static , using equations 1 and 2. 


E=- g r + h„ 


(i) 


Equation 1 is solved using gravity, g, and the estimated 
aircraft altitude, h ej( , from equation 2. 


h„i = 


I Pstatic \ 01903 

v 101325/ 


44333.7 


( 2 ) 


The total energy calculations are shown in blocks 220 and 
230. Total energy is filtered with filter 1 before it is differen- 
tiated with filter 2. Filter 1 is a 2" d order filter with a band- 
width of 1 .2 rad/s that is used to remove noise and measure- 
ment resolution errors from the estimated total energy before 
the signal is differentiated. Filter 2 is the combination of a 2 nd 
order filter and a differentiator. Filter 2 both filters and calcu- 
lates the rate of change of the input signal. 

Eneigy rate, E, is corrected to account for the energy added 
to the aircraft by the motor 210. The correction tenn was 
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found by differentiating equation [1] and solving for motor 
terms only. The resulting relationship is given in equation [3] . 


(3) 5 


6 

and the latitude and longitude of the chosen origin location, 
Lat^ and Lon,,. The origin location can be chosen to be the 
location of the UAV ground station. 

A history of previous aircraft positions and energy rates 
was collected in a first in first out queue, block 310, given by 
equation 9. 


Equation 3 is solved using an approximation of V motoJ . 
given in equation [4] 10 


^ motor — . . 


(4) 


15 

where T is the axial force from the motor and M is the 
aircraft mass. Thrust is calculated from the simple engine 
model given in equation [5] 


20 


0.64(l/irof,_i - throt Ir i m )T max 
s 2 + 0.16i + 0,64 


(51 


where throt^ is the 1 -second delayed throttle command 2 s 
and thiot trim is the trim throttle needed for level flight. The 
term 7 max is the maximum thrust generated by the motor. 
Equation [5] is essentially a gain on the throttle command 
delta from trim with a one-second time delay and a 2" d order 
filter with a 0.8 rad/s cutoff frequency. Values for T„ wt and 30 
throt . were initially taken from ground measurements but 
were later updated using flight data. 

Preferably, a limiter 250 is employed in order to ensure that 
the data calculated falls within certain boundaries, to assure 
that the data makes “real world” sense. 35 

Filter 3, 240, is used to calculate energy acceleration, E. 
Energy acceleration is normalized to produce normalized 
energy acceleration, E,„ using equation 6. 

40 

E ( 6 ) 

E n = — r 

ma x(£, 1 .0] 


Normalized energy acceleration is used for soaring control 45 
because normalization reduces the sensitivity of the control- 
ler to variations in thermal strength. 

The three filters used in the total energy calculations 
remove noise from the aircraft energy acceleration estimate 
but with a cost of approximately 2.5 seconds of time delay. 50 
This delay reduced the effectiveness of the E feedback chan- 
nel in the controller and is the primary reason for the devel- 
opment of a thermal state estimator described in blocks 20 
and 30 of FIG. 1. 

FIG. 3 depicts a more detailed explanation of the thermal 55 
identification block 20 of FIG. 1. 

An estimation of the thermal location, size, and strength is 
made in order to provide the controller with a faster indication 
of the aircraft state with respect to the thermal. The aircraft 
position, P, is calculated in block 300 from the latitude and 60 
longitude using equations 7 and 8. 


Px= (Lat- Lat 0 ) a 

(7) 

Py=(Lon-Lon 0 )a cos(Lato) 

(8) 


Equations 7 and 8 are solved using the equatorial radius of 
the Earth, a, the aircraft latitude and longitude, Lat and Lon, 


Px k Pr k L 

p h-i p h - 1 K-i 


Px k f/j 1 1 Px k-Lq + X P k-L; ■ 


(9) 


Entries to the queue given in equation 9 are made once per 
second, where k is the current measurement index. The length 
of the queue, L q , used in equation 9 was chosen to be 45. 

Estimation of the movement or drift of the thermal, block 
320, is important for achieving a good estimate of thermal 
position and radius. Drift is primarily caused by wind but 
early flight tests showed that thermals do not always drift with 
the same velocity or direction as the prevailing winds. Ther- 
mal drift is calculated by comparing the position of the ther- 
mal given by the first 20 entries of the queue with the position 
of the thermal given by the last 20 entries of the queue. To 
begin, the energy rate entries of the queue are shifted using 
equation 10 to prevent zero or negative values. 


9s = q\ 


1 0 
0 1 


0 0 


0 

0 

1 - min(#i ( : , 3)) + 0.01 


(10) 


Equation 1 0 produces the shifted queue, q^, where the term 
qj(:,3) represents the V d column of q r 

The thermal position given by the first or newest 20 rows of 
q,, is calculated using equation 1 1 . 


EupDi ~ 


20 

L q s (i, 1:2 )■?,(/, 3) 


20 

£ 9,(1, 3) 

;=i 


(ID 


The thermal position given by the last 20 rows of q, is 
calculated using equation 12. 


'upD2 ' 


N i 

L q,(l, 1 : 2 )q s (i, 3) 
<=«,— I 9 

N 

£ 9,(1, 3) 

;=/v,-i9 


( 12 ) 


The estimated drift velocity of the thermal, D, is calculated 
from equation 13 


EupD2 EupD] 

N„- 20 


(13) 
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The drift velocity is limited to +/-10 m/s and rate was 
limited to +/-0. 1 m/s 2 to prevent measurement errors in q, 
from causing large changes in D. The estimated drift velocity 
is then used to correct q, for drift to form a corrected queue, q^, 
using equations 14-15. 

q c (i,l:2)=(i-l)D+q s (i,l-2) (14) 

q c (i,3)=q„(i,3) (15) 

FIG. 4 shows the effect of the drift velocity in moving the 
positions in q^ to the positions given in q^. Typical data from 
soaring flight is used in this illustration. The result of the drift 
correction is to shift all the positions in the queue together to 
appear as though the thermal did not drift. The resulting 
queue, q^, can then be used to give a better estimate of the 
current thermal position and radius. 

Again, referring to FIG. 3, the estimated thermal position, 
P , is calculated, block 330, using the corrected queue, q c , as 
given in equation 16. 


Nq (16) 

^ q c (i, 1 : 2 )q c (i, 3) 2 

p = — 

“ P N q 

Z ?c(i. 3) 2 

i=l 


8 

Equation 19 is first solved using an initial value of 45 
meters for the estimated thermal radius, r , and then esti- 
mated thermal radius thereafter. The shape of the thermal 
velocity distribution is plotted in FIG. 5. 

5 A weighted error, err w , is calculated by taking the differ- 
ence between the predicted thermal velocity, w d , and the 
energy rate of the aircraft given in the 3 column of the queue, 
q c . The weighted error vector is calculated using equation 20. 

eT7 yXfy = ( W preJ})~qc(i’^))H{i) ( 20 ) 

The influence vector, H, is used to weight the more recent 
values ofq^ and w p „,^ higher than the older values. Values for 
H are plotted in FIG. 6. 

The sum squared error, err su , is calculated using equa- 
tion 21. 


err w -err w n\\ 

err sumsq = 77 

/V <? 

20 

Equation 2 1 gives a measure of the fit between the assumed 
thermal shape given in equation 1 9 and the measured thermal 
velocity stored in the queue. The sum squared error for a 
perturbed radius, err sumsq , is calculated in a similar way using 
25 equations 22 through 24. 


Equation 1 6 differs from equations 1 1 and 1 2 in that the E 
terms of the queue are squared. This makes the resulting 
estimated thermal position move closer to areas of strong lift. 
Equation 1 6 represents the centroid of the aircraft energy rate 
during the previous N ? seconds. This method for estimating 
thermal position is easy to implement, does not require tun- 
ing, can be used for any thermal size, and does not have a high 
computational cost. Drawbacks associated with this method 
are that it is sensitive to the length of the queue and it usually 
biases the estimated thermal position toward the center of the 
measurement set. 

Thermal velocity, w up , is approximated, block 340, from 
the aircraft energy rate history found in q c . The equation used 
to estimate w up is given in equation 17. 

w^=l.lmax(? c (:,3)) (17) 

Equation 1 7 uses the assumption that the maximum ther- 
mal velocity is 10% greater than the maximum climb rate 
experienced by the vehicle during the previous N seconds. 
The thermal velocity is rate limited to not increase faster than 
0.025 m/s 2 and not to decrease faster than 0.0 1 5 m/s 2 as a way 
to retain old information about the thermal strength that has 
been passed from the queue. The estimated thermal strength 
is used in addition to q c to estimate the thermal radius. 

The radius of the thermal, block 350, is estimated using a 
gradient descent method to fit an assumed thermal shape to 
the data found in the queue. The estimation process starts with 
the calculation of the distance between each position in q c to 
the estimated thermal center, P using equation 18. 

S(i)=,|( ?c (U)-i ) „ p (i) 2 +(? c ft2)-P„ p (2)) 2 (18) 

The vector, S, is used to calculate the predicted thermal 
velocity of each point in q c using the thermal velocity distri- 
bution shape given in equation 19. 


f « ) 2 (19) 

= (»v + 0.5)e _ o.5 


(j T (22) 

Wpred = (Wup + 0.5)e -0.5 


erf^{u- pretr q c (:,3))H (23) 


. err w • err w (24) 

err sumsq — ^ 

The perturbed radius, r , used in equation 22 is defined by 
equation 25. 

'V=v+6 r (25) 

The stepsize, 8, is set to 0.5 m for this illustration. 

The sum squared error for the current thermal radius and 
the sum squared error for the perturbed thermal radius are 
used to calculate a simple gradient as given in equation 26. 


e ^^sumsq err sumsq (— 6 ) 


The gradient is then used to calculate the thermal radius for 
the next frame using equation 27. 

'Wi =L *r+ r „ Pt (27) 

The learning rate, L, is set to 10 for this illustration. 

Circle guidance block 30 of FIG. 1 is described in more 
detail below. Guidance calculations are used to determine the 
steady state turn rate needed to fly in the estimated thermal 
and to produce position and velocity errors to the soaring 
controller. The aircraft is commanded to fly a circular flight 
path having a radius, r cmd , determined by equation 28. 

r„, Itr 0.65r„ p (28) 

The scale factor 0.65 is chosen because it defines a flight 
path radius that is usually small enough to be within the 


30 


35 


40 


45 


50 


55 


60 


65 
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thermal “core” but laige enough to avoid the high aircraft sink 
rate associated with high bank angle turns. The steady state 
turn rate, needed to fly at an airspeed ofV along a circular 

path given by r cmd is given in equation 29. 


-180V (29) 

= — ■ — 
cmd 

The negative sign found in equation 29 is used to ensure 
that the aircraft always circles to the left while soaring. This is 
done for simplicity and to give the aircraft ground personnel 
a quick visual indication of the aircraft mode during flight 
tests. For this embodiment of the invention, flight plans were 
designed to command right turns when the aircraft is search- 
ing for thermals and left turns when the aircraft is in soaring 
mode. Position error, err p , is generated using equation 30. 

err p =r cmcr S( 1) (30) 

Position error is differentiated and filtered to create veloc- 
ity error, err v using the transfer function given in equation 3 1 . 

1.69s (31) 


The guidance calculations given inequations 28-31 calcu- 
late the position error, velocity error, and turn rate command 
information that is used by the soaring controller to maintain 
a circular flight path that was centered on the thermal. 

Referring to FIG. 7, the controller block 40 from FIG. 1 is 
described in further detail. The soaring controller 700 archi- 
tecture is loosely based on a thermal centering method given 
for glider pilots in “Cross Country Soaring” by Reichmann, et 
al., Soaring Society of America, Inc., Hobbs, N. Mex., ISBN 
10883813-01-8, 1993, which is incorporated herein by refer- 
ence. The method states three rules; 

1. As climb improves, flatten the circle (approx. 15-20°) 

2. As climb deteriorates, steepen the circle (approx. 50°) 

3. If climb remains constant, keep constant bank (approx 
25-30°) 

These rules are used as a guideline for the soaring control- 
ler design. Energy acceleration feedback is used to implement 
rules 1 and 2 because the rules given by Reichmann, et al. are 
to be applied to the “changes in rate of climb.” Rule 3 is 
implemented with steady-state turn rate. The soaring control- 
ler 700 inputs consist of normalized energy acceleration, E„, 
position error, err^, velocity error, err v , and steady-state turn 
rate, The soaring controller 700 output consists of a turn 
rate command, that is sent to the autopilot. The con- 
troller gains, k^kj, and k 3 , 710, 720, and 730 respectively, are 
tuned by first setting k 2 720 and k 3 730 equal to zero. The 
energy acceleration gain k x 710 is tuned to produce a simu- 
lated aircraft response that followed the rules outlined by 
Reichmann, et al. when subjected to a variety of thermal 
strengths and sizes. The thermal model given in the reference 
is implemented into a hardware in the loop simulation for 
controller 700 evaluation and gain tuning. The k 2 720 and k 3 
730 gains are timed to make the path of the aircraft more 
circular and to improve the damping of the controller 700 
during soaring flight in a variety of thermal sizes and 
strengths. The resulting control gains were chosen to be 50, 
0.4, and 0.165 for k, , k 2 , and k 3 respectively. Limiters 740, 
750, 760, and 770 are used to limit the inputs into the con- 
troller 700 to ensure that one input does not dominate the 
aircraft’s control. The output turn rate command is also lim- 
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ited by a turn rate limiter 780 to ensure the aircraft doesn’t 
turn too abruptly and a limiter 790 to ensure the aircraft only 
turns left during soaring. Alternatively, the mm rate command 
output of the soaring controller can be used as a bank angle 
5 command if desired. 

FIG. 8 shows a mode logic diagram that depicts the mode 
switching block 50 from FIG. 1 in greater detail. 

In general, mode logic is used to determine when to switch 
between searching flight and soaring flight. The mode logic 
to described in this section is significant because it determines 
how a thermal is detected. The mode logic is broken into three 
components. Thermal detection logic 800 mode switching 
uses aircraft response to detect thermals upon encountering 
thennals. An encounter with a thermal by the aircraft pro- 
15 duces large changes in energy rate and energy acceleration 
because of the atmospheric eneigy imputed to the aircraft as 
it traverses the thermal. During flight through a thermal, the 
energy acceleration changes sign and becomes negative when 
the energy rate peaks. The thermal detection logic is tuned to 
20 trigger when the thermal is imputing significant eneigy to the 
aircraft and the peak value lias been reached. Thus, the ther- 
mal detection logic 800 provides a basic crest determination, 
above a selected threshold, for the thermal. 

The thermal rejection logic 810 provides a decision to 
25 switch from soaring flight to searching flight using energy 
rate. Thermal rejection 81 0 is needed to switch from soaring 
flight to searching flight when the thermal is no longer pro- 
viding sufficient energy to the aircraft. Two energy rate based 
criteria are used by the thermal rejection logic 810. First, 
30 thermals are rejected if the smoothed output of filter 4 is less 
than -0.5 m/s. This criteria is met if the aircraft encountered 
light sink for a period of approximately 5 seconds. The second 
criteria, based on the output of filter 5, is triggered if the 
aircraft has not experienced sufficient energy rate for approxi- 
35 mately 10 seconds. Additional criteria, not shown, may be 
used to disengage from soaring mode. The additional criteria, 
for example, are upper altitude limit exceeded, lower altitude 
limit exceeded, lost link, and soaring mode disabled by 
ground station operator. 

40 The third component is the latch mechanism 820. The latch 
mechanism 820 receives inputs from the thermal detection 
logic 800 and the thermal rejection logic 810 and provides the 
actual “switch” command from soaring flight to searching 
flight and visa-versa to the aircraft. 

45 What is described are specific examples of many possible 
variations on the same invention and are not intended in a 
limiting sense. The claimed invention can be practiced using 
other variations not specifically described above. 

What is claimed is: 

50 1. A method for locating and using thennals to improve 

endurance of an uninhabited air vehicle comprising: 

calculating eneigy rate and energy acceleration from total 
energy of the air vehicle using static pressure around the 
air vehicle, tme air speed of the vehicle, and air vehicle 
55 engine input; 

calculating a thermal radius, velocity, and position using 
the energy rate and a location of the air vehicle; 
calculating air vehicle position error and velocity error 
using location of the air vehicle, the thermal radius, the 
60 thermal velocity, and the thermal position; 

calculating a soaring turn rate for the air vehicle using the 
position error, the velocity error, and a steady-state turn 
rate required to fly within the thermal when the air 
vehicle meets a selected climb rate; 

65 calculating a final turn rate for the air vehicle through mode 

switching using the energy acceleration, the soaring turn 
rate, and a waypoint tracking turn rate; and, 
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controlling the air vehicle to guide the air vehicle to the 
thermal to stay within the thermal when the air vehicle 
meets a selected climb rate. 

2. The method of claim 1, wherein the calculating energy 
rate and energy acceleration step further comprises: 

calculating an engine input effect on energy rate using the 
true air speed and air vehicle engine input; 
calculating a total energy using air vehicle velocity and 
static pressure around the air vehicle; 
filtering the total energy to remove noise and measurement 
resolution errors; 

differentiating the filtered total energy to obtain a rate of 
change of input; 

calculating energy rate by adding the engine input effect to 
the rate of change of input; and, 
calculating the energy acceleration by differentiating the 
energy rate. 

3. The method of claim 1, wherein calculating the thermal 
radius further comprises: 

calculating the thermal velocity by building a history of 
previous aircraft locations using previous aircraft posi- 
tions and energy rates to obtain a thermal velocity esti- 
mation; 

calculating thermal drift correction using the history of 
previous aircraft locations and energy rates; 
calculating thermal location by calculating a center of the 
thermal position using an aircraft location obtained after 
thermal drift correction; and, 
calculating thermal radius using the location of the air 
vehicle, the thermal velocity calculation, the thermal 
drift correction, the aircraft location after thermal drift 
correction, and the center of the thermal position. 
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4. The method of claim 1, wherein the mode switching 
comprises: 

using normalized energy acceleration and energy rate to 
make a determination of when the air vehicle reaches a 
5 crest of the thermal; 

using energy rate to make a lift determination of when lift 
caused by the thermal falls below a selected level; and, 
calculating a thermal stopping point wherein the air vehicle 
leaves the thermal using the crest determination and lift 
to determination. 

5. The method of claim 1, wherein the controller step 
comprises controlling an air vehicle turn rate and velocity to 
guide the air vehicle to an approximate center of the thermal. 

6. The method of claim 5, wherein controlling air vehicle 
15 turn rate comprises the steps of: 

flattening the turn rate as an air vehicle climb rate 
increases; 

steepening the mm rate as the air vehicle climb rate 
decreases; and, 

20 maintaining a constant turn rate as the air vehicle climb rate 

is constant. 

7. The method of claim 6, wherein the flattening and steep- 
ening steps use energy acceleration gain data. 

8. The method of claim 5, wherein the controlling step 
25 comprises inputs of normalized energy acceleration, position 

error, velocity error, and steady-state mm rate. 

9. The method of claim 8, wherein the inputs are limited to 
prevent a single input from dominating the result. 

10. The method of claim 9, wherein the controlling step 
30 only allows the air vehicle to turn left. 



